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Abstract. We present the analysis of broad absorption lines (BALs) seen in the spectrum of the z cm ~4.92 QSO 
SDSS J160501. 21-011220.0. Our high spectral resolution UVES spectrum shows two well detached absorption line 
systems at z a b s = 4.685 and 4.855. The system at z a ^ s — 4.855 covers the background source completely suggesting 
that the gas is located outside the broad emission line region. On the contrary the system at z a bs= 4.685, which 
covers only on the continuum source, has a covering factor of the order of 0.9. Physical conditions are investigated 
in the BAL system at z a b s = 4.855 using detailed photoionization models. The observed H I absorption line together 
with the limits on C II and Si II absorptions suggest that 16 < log JV(H i) (cm -2 ) < 17 in this system. Comparison 
with models show that the observed column densities of N v , Si IV and C IV in this system require that nitrogen 
is underabundant by more than a factor 3 compared to silicon if the ionizing radiation is similar to a typical QSO 
spectrum. This is contrary to what is usually derived for the emission line gas in QSOs. We show that the relative 
suppression in the N v column density can be explained for Solar abundance ratios or abundance ratios typical of 
Starburst abundances if an ionizing spectrum devoid of X-rays is used instead. Thus, if the composition of BAL 
is like that of the emission line regions it is most likely that the cloud sees a spectrum devoid of X-rays similar 
to what we observe from this QSO. This is consistent with the fact that none of our models have high Compton 
optical depth to remove X-rays from the QSO. Similar arguments lead to the conclusion that the system at z a t, s = 
4.685 as well is not Compton thick. Using simple Eddington arguments we show that the mass of the central black 
hole is ~ 8 x 10 8 M@. This suggests that the accretion onto a seed black hole must have started as early as z = 
11. 

Key words. Quasars: absorption lines - Quasars: individual: SDSS J160501. 21-011220.0 



1. Introduction 

Broad absorption line systems (BALs) seen in the spec- 
tra of QSOs are characterized by absorption features 
with large velocity widths ( Av ~ a few 1000 km s -1 ) 
and, usually, high ionization states (Turnshek ct al. 1988; 
Weymann et al. 1991). The gas giving rise to BALs is be- 
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lieved to be material ejected by the quasar but still located 
very close to the central regions. Thus studying metallic- 
ities in the absorbing gas is a direct probe of the chem- 
ical enrichment in the very central regions of AGNs (see 
Hamann & Ferland 1999 for the review). The associated 
absorption line systems by and large show large metallici- 
ties with a tendency of nitrogen being over-abundant with 
respect to oxygen and other a— process elements (e.g. 
Petitjean et al. 1994; Korista et al. 1996; Hamann 1997; 
Petitjean & Srianand 1999). Detailed studies of emission 
line properties of z > 4 QSOs have revealed Super-Solar 
metallicities in the line emitting gas in broad line regions 
(BLR), suggesting very rapid star formation process (e.g. 



2 



Neeraj et al: BAL QSO SDSS J160501.21-011220.0 



3 




2 6600 6800 7000 7200 7400 

Sh 

X 




7600 8000 8400 8800 9200 

Observed Wavelength (A) 

Fig. 1. UVES spectrum of the z om = 4.92 QSO SDSS J160501. 21-011220. Two well detached BAh features are seen 
covering the redshift ranges 4.62—4.75 and 4.83—4.88. We denote them as blue and red systems respectively. The 
expected positions of different lines (with sub-scripts r and b denoting red and blue respectively) are marked with 
solid lines. The horizontal dotted line in the top panel gives the mean transmitted flux computed over the Lya forest 
in the wavelength range 6550—6800 A that is not affected by the BAL absorption lines. 



Hamann & Ferland 1992; Dietrich et al. 2003). However, 
no such metallicity estimates at such high redshifts based 
on absorption lines using echelle spectra are available. 

Unlike in the case of most normal galaxies, it is believed 
that the high metal enrichment usually observed in QSOs 
has taken place in less than a few 10 8 years through rapid 
star formation similar to what is expected to happen in the 
center of ellipticals (Hamann & Ferland 1993; Matteucci 
& Padovani 1993). Thus, studying BAL QSOs at z cm > 5, 
where the age of the Universe is t< 1.07 x 10 9 yr (for 
H = 75 km s" 1 Mpc" 1 , fl A = 0.70, Ct m = 0.30), is im- 



portant for understanding the star-formation history at a 
time close to the epoch of reionization and possibly also 
for constraining cosmological parameters. In the initial list 
of SDSS QSOs (Fan et al. 2000) there are few QSOs at 
z > 4 with associated absorption lines. We obtained a 
UVES spectrum of the QSO SDSS J160501. 21-011220.0, 
the highest redshift BAL QSO known at that time, with 
the aim to study the outflowing gas in detail. 

In order to derive a realistic estimate of the absorbing gas 
metallicity one needs to have a good handle on the ion- 
ization corrections. Associated systems being close to the 
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QSO are most probably ionized by the QSO light rather 
than the diffuse intergalactic background. However the 
ionizing spectrum from the BAL QSOs is poorly known. 
It is known that BAL QSOs are under-luminous in X-rays 
(Bregman 1984; Singh et al. 1987; Green & Mathur 1996). 
Recent, deep Chandra observations have shown that the 
optical to X-ray spectral index, aox-, measured for BAL 
QSOs is systematically lower than that of non-BAL QSOs 
(Green et al. 2001). It is not clear whether this X-ray weak- 
ness is due to the QSO being intrinsically X-ray quiet (i.e., 
this is a property of the central engine) or due to line-of- 
sight absorption (i.e., radiative transfer effects). Using a 
complete sample of optically selected QSOs, Brandt et al. 
(2000) have shown a significant anti-correlation between 
aox and the total equivalent width of the C IV absorption 
lines. In addition it is noticed that the emission line prop- 
erties of the BAL and non-BAL QSOs are very similar 
even though the observed X-ray properties differ signifi- 
cantly (Weymann et al. 1991; Korista et al. 1993). Thus, 
it is most likely that the relative X-ray weakness of BAL 
QSOs is due to intrinsic absorption. 

It must be noted that all existing models of BAL sys- 
tems assume standard QSO spectrum. However estimated 
metallicities will be different if the intrinsic aox is differ- 
ent from that of a typical QSO spectrum. Indeed, Srianand 
& Petitjean (2000) have shown that the inferred metallic- 
ities decrease with decreasing value of aox- Also, if the 
absorption is responsible for the suppression of X-rays 
then it is important to investigate whether it is caused 
by the gas responsible for the BAL troughs (common ab- 
sorbers) or by a distinct gas component. It is likely that 
the common absorber picture is correct for the BAL to- 
wards PHL 5200 (Mathur et al. 1995) and that a distinct 
Compton thick screen is required for the BAL towards 
PG 0946+301 (Gallagher et al. 1999; Mathur et al. 2000; 
see also Kraemer et al. 2002). 

In this paper we present a detailed analysis of two 
well detached BAL outflow components that are seen in 
the spectrum of QSO J160501. 21-011220.0 (z em = 4.92). 
Details of the observations are discussed in Sect. 2 which 
is followed by the description of the BALs and column 
density estimates in Sect. 3. In Sect. 4 we investigate the 
physical conditions of the gas. The results are summarised 
in Sect. 5. 

2. Observations and data reduction 

The Ultraviolet and Visible Echelle Spectrograph (UVES; 
Dekker et al. 2000), installed at the ESO VLT 8.2 m tele- 
scope, unit Kueyen, on Mount Paranal in Chile, was used 
on June 15-17, 2001, to obtain high resolution spectra of 
QSO J 160501.21-011220.0. A non-standard setting with 
cross-disperser #4 and central wavelength 8420 A was 
used in the Red arm of UVES. Full wavelength coverage 
was obtained this way from 6514 to 8300 A and from 8494 
to 9300 A accounting for the gap between the two Red- 
arm CCDs. The CCD pixels were binned 2x3 (namely, 
twice in the spatial direction and three times in the disper- 



sion direction) and the slit width was either fixed to 1'.'2 
or l'/5, yielding an overall spectral resolution R ~ 30000. 
The total integration time 3h30min was split into 3 ex- 
posures. The data were reduced in the dedicated context 
of MIDAS 1 , the ESO data reduction system, using the 
UVES pipeline (Ballester et al. 2000) in an interactive 
mode. The main characteristics of the pipeline is to per- 
form a precise inter-order background subtraction for sci- 
ence frames and master flat-fields, and an optimal extrac- 
tion of the object signal rejecting a number of cosmic ray 
impacts and subtracting the sky spectrum simultaneously. 
The pipeline products are checked step by step. The wave- 
length scale of the spectra reduced by the pipeline was 
then converted to vacuum-heliocentric values and individ- 
ual 1-D exposures scaled, weighted and combined alto- 
gether using the NOAO onedspec package of the IRAF 2 
software. During this process, the spectra were rebinned 
to 0.08 A pix -1 . In order to increase the signal-to-noise 
ratio in this near-IR part of the optical range (and as 
the lines of interest are broad), we applied a Gaussian fil- 
ter smoothing with two pixel FWHM and get an effective 
spectral resolution of ~ 12 km s _1 . The resulting signal- 
to-noise ratio per pixel is of the order of 20 or more over 
most of the wavelength range considered here. The final 
spectrum is shown in Fig. ^ 

The continuum fitting in the red side of the Lya emis- 
sion has been performed using a smooth low order poly- 
nomial considering only the absorption free regions. In the 
Lyman-a forest, we approximated the QSO spectrum by a 
powerlaw, f v oc v~ a e , that we fitted on the low-dispersion 
SDSS spectrum made available to us by Dr. Fan. The nor- 
malization of the power-law was done using the flux in the 
region between the C iv and Si iv emission lines. In our 
normalised spectrum the mean Lya transmission flux in 
the redshift range z=4.39 to z=4.60 (that is not contam- 
inated by any BAL absorption lines) is 0.15. This is con- 
sistent with the values derived in the same redshift range 
using non-BAL QSOs (Becker et al. 2001). 

3. Broad absorption line systems 

Based on the unabsorbed C IV emission line we estimate 
the emission redshift of the QSO to be z cm =4.92. Two 
well detached absorption line systems are identified based 
on the C iv and Si iv absorptions over the redshift ranges, 
z abs =4.83-4.88 and z abs = 4.62-4.75, respectively. For 
mean redshifts z a b s = 4.855 and 4.685, this corresponds to 
ejection velocities of «3,330 km s _1 and ps 12,160 km s -1 . 
The spectrum in the Lyman-a forest is heavily affected by 
the high Lya opacity of the intergalactic medium. Even 
though excess absorption is seen at the expected positions 

1 MIDAS: Munich Image Data Analysis System, trademark 
of the European Southern Observatory (ESO) 

2 IRAF: The Image Reduction and Analysis Facility is dis- 
tributed by the National Optical Astronomy Observatories, 
which is operated by the Association of Universities for 
Research in astronomy, Inc. (AURA), under cooperative agree- 
ment with the National Science Foundation. 
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Relative velocity (kms 1 ) 

Fig. 2. Velocity plots of the red absorption line system 
centered at z a b s = 4.855. We identify 5 distinct compo- 
nents based on the Si iv and C iv profiles (marked with 
vertical dotted lines). Detection of Lya line in this system 
is unambiguous as the profile just follows the well detached 
Si iv absorption line at least from to 600 km s . 

of BAL absorption lines (see Fig. ^) it is difficult to es- 
timate column densities in this region. Therefore, in the 
following sections we use only lines observed on the red 
side of the Lya emission line to investigate the physical 
conditions in the outflowing gas. 

3.1. The red component at z a h s = 4.855 

The absorption trough in this component is ~ 2000 
km s _1 wide. The absorption profiles of different tran- 
sitions centered at z a hs = 4.855 are shown in Fig. [2J The 
absorption profile of the Si iv doublet that is compara- 
tively less saturated than that of the C iv doublet suggests 
the presence of 2 relatively narrow components denoted 
by 11 and 12 and 3 broad subcomponents denoted by a, b 
and c. While Lya and C iv are detected in all these com- 
ponents, N v absorption is detected only in the broad 
subcomponents (see Fig.[2jl. 

All the obviously saturated absorption lines have zero 
residual (C IV and Lyman-a, see Fig.|2J). The covering fac- 
tor of the gas is therefore unity. In case of complete cover- 
age the observed residual intensities Ri and R 2 in the first 
and second member of the doublet are related as Ri^R?, 
(see Srianand and Shankarnaraynan 1999). Indeed, even 
though part of the N v absorption line is contaminated by 
atmospheric absorption (see Fig.^J, the velocity ranges in 
the profiles that are free of blending are consistent with 
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Fig. 3. Velocity plots of N v doublets centered at z a b s = 
4.855. Ri and R2 are the residual intensities in the first 
and second lines of the doublet. Vertical dotted lines show 
the subcomponents 11, a, b, c and 12 respectively. In the 
lower two panels the dotted profiles give the atmospheric 
lines seen toward QSO 1122—1628. In the top panel we 
overplot the two N v profiles after appropriately scaling 
(see text) the N vA1242 assuming complete coverage. The 
good matching of profiles suggests that N v absorption 
line covers the background source completely. 

complete coverage (see top panel in Fig.^J. As emission in 
this wavelength range is contributed by both continuum 
and broad emission lines the corresponding absorbing gas 
has to be located outside the broad line region (BLR). As 
the N v lines are weak we can derive the N v column densi- 
ties. However, due to saturation effects, we can determine 
only lower limits for C IV and H 1 column densities. It 
can be seen on Fig. [21 that Si IVA1402 is only partly satu- 
rated which means that, although very uncertain, we can 
estimate the corresponding column density. We use also 
the Si IVA1402 profile as a template to determine upper 
limits on the column densities of non detected species (see 
below) . 

The column density in each velocity pixel was obtained 
using the relation 

N(v) = 3.768 x 10+ 14 r/fX cm" 2 km _1 s (1) 

where, r, / and A are the optical depth, oscillator strength 
and rest wavelength respectively. The total column den- 
sity is obtained by integrating N(v) over the velocity 
range —1120 to 1160 km s . The results are presented 
in Tabled Rest wavelengths and oscillator strengths used 
here are taken from Verner et al. (1994). The two lines 
of the C iv doublet are blended and we use an effective 
oscillator strength, / = f± + /a, and mean wavelength, 
(Ai + X%)/2, to compute the lower limit on the column den- 
sity. In the case of N v, we use the N vA1238 line for com- 
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Fig. 4. Velocity plots for the blue component centered at 
z a bs — 4.685. Second member of the doublets is overlaid as 
a dotted line. Vertical dashed lines mark the velocity range 
over which Si iv absorption line is seen. Note that C IV ab- 
sorption is seen over a larger velocity range. The non-zero 
residual flux at the expected position of Lycn absorption 
is used for obtaining the upper limit on the H I column 
density. 

ponent a and N vA1242 line for components b and c that 
are free of blending. Contamination by atmospheric fea- 
tures (see Fig. is corrected by carefully masking the at- 
mospheric contamination using the normalised spectrum 
of Q 1122-1628. 

Singly ionized species such as Si II, C n and Al n are ab- 
sent and the spectrum at the expected position of Al III 
lines is unfortunately very noisy. As said above, we evalu- 
ate upper limits for column densities of the latter species 
using the Si IVA1402 profile (which is not completely sat- 
urated and least affected by atmospheric contamination) 
extending from —1120 to 1160 km s _1 as a template. 
The scaling factor k — [N fX]x+ /[N fX\template between 
the two optical depths Ttempiate and t x + for species X + 
is then obtained by minimizing 

a = S (r x + —k* Ttempiate) 2 (2) 

For Si II and C II, k comes out to be 0.031 and 0.035 
respectively. The corresponding upper limits on the col- 
umn densities are given in Table [TJ Note that the error 
in the column density of N v is mainly due to continuum 
placement uncertainties. We note that most of the atmo- 
spheric absorption seen in the expected wavelength range 
have consistent equivalent width. This suggests that our 
continuum fitting does not underpredict the absorption in 
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Fig. 5. Region of the spectrum containing the Si iv trough 
of the red BAL component at z a b s ~ 4.855. Note that 
systems 11 and 12 have a velocity separation equal to that 
of the Si IV doublet splitting suggesting the presence of 
line-locking in the system. 

Table 1. Column densities estimated using apparent op- 
tical depth method with complete coverage. 





species 


logN (cm -2 ) 


4.855 


H i 


>15.67 




C n 


<14.21 




C iv 


>15.64 




Si ii 


<14.13 




Si iv 


15.33: 




N v 


14.63-14.71 


4.685 


H i 


< 16.00 




C iv 


>15.94 




Si ii 


<14.73 




Si iv 


>15.40 



Table 2. Chemical compositions considered in the models: 





Solar 




Starburst* 




Metal 


1.0 z 


1.0 Z SB 


0.5 Z SB 


0.20 Z ss 


C 


-3.44 


-3.90 


-4.22 


-4.63 


N 


-4.03 


-4.60 


-5.34 


-6.51 


O 


-3.13 


-3.00 


-3.30 


-3.70 


Si 


-4.44 


-4.43 


-4.74 


-5.14 



* This form of Starburst abundances correspond to chemical 
evolution model M5a of Hamann & Ferland (1993). 



the N V region. No flux is detected at the expected posi- 
tion of Si IHA1206 (see Fig.^l. However, we could not use 
this to derive a lower limit on Si in as the same region is 
contaminated by a probable N v absorption line from the 
blue component (see Fig.0J|. 

Finally we notice that the velocity separation between the 
two narrow components 11 and 12 of the red component is 
very close to the Si iv doublet splitting. This is shown in 
Fig. [SJ This is consistent with growing evidence for line- 
locked flows and radiative acceleration in BAL systems 
(see Srianand et al. 2002 and references therein). 
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3.2. The blue component at z a b s = 4.685 

The absorption profiles of different transitions centered 
at z a bs = 4.685 are shown in Fig. Doublet partners for 
C iv and Si iv are partially blended together. The non-zero 
residual flux for the blended C ivAA1548,1550 corresponds 
to a covering factor, f c >0.90. It is interesting to note that 
the absorption is well detached from the emission line pro- 
file. Thus the partial coverage reflects either ~ 10% con- 
tribution from the scattered light or partial coverage of 
the continuum source. Lower limits to column densities 
for the Si ivAA1393,1402 and C ivAA1548,1550 were ob- 
tained assuming full coverage (see Table P). Doing other- 
wise would only increase these limits. Note that C iv ab- 
sorption line extends over a larger velocity range as com- 
pared to Si IV absorption profile (see Fig.0J). Like in the 
red component, singly ionized species Si II and C II are 
not detected which indicates that the neutral hydrogen 
column density cannot be large (remember that metallic- 
ity in this kind of gas is usually large) . N v may be present 
but its redshifted position in the Lya forest overlaps with 
that of H I and Si ill from the red component (see Fig.QJ. 
The mean transmitted flux at the N v position is consis- 
tent with saturated absorption line. 

In the Fig. 0] the intervening Lyman-a forest is clearly 
observed over the whole range corresponding to the BAL 
H I Lyman-a trough. This means that the H I optical 
depth cannot be large except if the covering factor of H I is 
much smaller than that of C iv. This is not impossible (see 
e.g. Srianand et al. 2002) but would be surprising given 
the non detection of both C n and Si II. By comparing the 
mean transmission in this region with that devoid of BAL 
absorption lines we estimate that t(v) for H I is less than 
2. If we use the velocity range covered by the C iv profile 
(~ 4000 km s" 1 ), we see that the H I column density can 
not much larger than 10 16 cm" 2 . 

4. Discussion 

4.1. Physical conditions in the outflowing gas 

From the discussion in the previous Section, it is clear that 
constraints are not strong enough to probe the physical 
conditions in the blue component (z aDS = 4.685). We there- 
fore mainly concentrate on the red component (z aD s = 
4.855). Recent Chandra observations (Vignali et al. 2001) 
fail to detect X-rays from the QSO SDSS J160501.21- 
011220.0. As discussed before, the X-ray non-detection 
could be either due to large X-ray absorbing column den- 
sity or a consequence of the QSO being intrinsically X-ray 
quiet. Indeed, there are already some reports that the X- 
ray continuum shapes of QSOs may evolve at z >2.5 (e.g. 
Vignali et al. 1999; Blair et al. 2000; Vignali et al. 2003). 
Assuming that the lack of X-rays in this source is due to 
absorption, Vignali et al. (2001) infer that the X-ray non- 
detection for this object implies a total hydrogen column 
density of N(H(total)) > 5.0 x 10 23 cm" 2 , characteristic of 
low-ionization BAL QSOs at low redshift. Note that the 
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Fig. 6. The spectral energy distribution (SED) of the 
ionizing radiation. The continuous and dashed curves cor- 
respond to the SED of the MF and BAL spectrum re- 
spectively. The vertical dotted lines give the ionization 
potential of the different species we are interested in this 
study. 

limit on the X-ray flux further implies a small optical to 
X-ray spectral index, a ox < —1.82, for this object. 
The absence of singly ionized species C n and Si II sug- 
gests that H I is optically thin at the Lyman limit. Thus, 
to probe the physical conditions in the red component, we 
run grids of photoionization models using Cloudy (Ferland 
1996) in the range log N (H i) (cm" 2 )~ 10 16 to 10 17 
cm" 2 , considering the gas being ionized by either an unat- 
tenuated QSO spectrum given by Mathews & Ferland 
(1987) (hereafter MF spectrum; see Fig. EJ) or a modified 
spectrum with little X-rays (hereafter BAL spectrum; see 
Fig. EJ) mimicking attenuation of a typical QSO spectrum 
by a large column density of ionized gas (see below) . The 
models are run for different chemical composition, either 
the Solar one or the so-called Starburst one (see Table |2J). 
In all the models, the calculations are stopped when the 
total neutral hydrogen column density reaches the limit 
we set (log N(H i) (cm" 2 ) = 16 or 17). 
The input parameters of the models, that reproduce the 
observed ratios as well as the column densities of in- 
dividual species, are given in Table The correspond- 
ing total hydrogen column density in the allowed range 
of ionization parameters is given in the last Column 
of the Table Results for some selected models with 
different ionizing spectra and chemical composition for 
log iV(H i) (cm" 2 )) = 17 are given in Fig.0 In this fig- 
ure we have plotted the column density of different species 
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allowed range in log U from log(N(H(total)))* 

log N(H i)* ratios observed N 



Model parameters 

Spectrum enrichment Metallicity(Z) 

MF Solar 0.5Z Q ,[Si/N]=3.O[Si/N] 17.00 

BAL " I.OZq 16.00 

1.5Z 

0.05© 17.00 

0.1Z Q 

0.2Z q 

I.OZq 

Starburst l.OZss 16.00 

5.0Z SS 
9.0Z SS 

0.2Z SS 17.00 

0.5Z SS 

l.OZss 



* in cm 2 

as a function of ionization parameter. The horizontal dot- 
ted lines mark limits on the observed column densities of 
Si II , Si iv , C II and C iv . The minimum and maximum 
of the allowed range for log U are marked with verticle 
dashed lines. The former is obtained using either of the 
ratios AT(Si n)/JV(Si iv) or N(C n)/N(C iv) and the lat- 
ter is obtained from the observed 2V(N v)/N(Si iv) ratio 
for metallicities assumed in the model. The range of ion- 
ization parameters for which the observed N v column 
density is reproduced is indicated as a verticle shaded re- 
gion. 

We note that, if the H I column density is close to the ob- 
served value (i.e. few 10 16 cm -2 ) one needs overall metal- 
licities larger than Solar to explain the observed column 
densities. We therefore mainly concentrate on models with 
A^(H i)= 10 17 cm~ 2 that will give a conservative lower 
limit on the metal enrichment in the system. 

If the gas is ionized by a MF type spectrum and has a 
Solar composition, it can be seen from panel (a) in Fig. [7| 
that the observed constraints on the ion ratios and the 
N v and C iv column densities can be consistently repro- 
duced for log U ~ —1.9 (see the shaded region). However 
the predicted column density of Si IV is lower than the 
observed value by more than an order of magnitude. It 
is clear that the observed value can be reproduced only if 
[Si/H]>[Si/H] . Increasing the overall abundance actually 
reduces the required overabundance of Si with respect to 
N. For example if one assumes Z = 0.5 Z Q the observed col- 
umn densities can be reproduced consistently with lower 
ionization parameter and a [Si/N] ratio higher than Solar 
by a factor of 3 only (see Table. Larger metallicities 
are not possible however as C n would then be detectable. 
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The enhancement of AT(Si iv) with respect to N(N v) 
can be achieved by either enhancing the silicon metallicity 
with respect to the nitrogen one as discussed above or by 
suppressing the N IV ionizing photons. The first possibility 
can be achieved in a chemical enrichment model where the 
secondary production of N has not yet begun to dominate. 
The latter possibility naturally arise if the cloud sees an 
attenuated X-ray spectrum, similar to what we observe 
directly for SDSS J160501. 21-011220.0. We explore this 
possibility in the following. 

We model the absorbed spectrum (called BAL spectrum 
in the following) as the composite of a black body spec- 
trum with temperature T = 150,000 K, plus a power-law 
spectrum with a uv — —0.5 and a x = —1.0 and a relative 
scaling a ox — —2.0. From Fig. El it is clear that for a given 
number of hydrogen ionizing photons this spectrum has 
less photons to ionize Si III and N iv compared to the MF 
spectrum. Thus a given ionization state will be produced 
by higher values of log U and hence higher total hydro- 
gen column densities. However the major difference is the 
presence of a soft X-ray excess in the MF spectrum that 
implies more N IV ionizing photons for a given number of 
Si in ionizing photons compared to the BAL spectrum. 
For a column density of log N (H i) (cm _2 )«17.00, ob- 
servations can be reproduced by models with metallicities 
0.05Z Q <Z<0.2Z Q (Tabled. This is also demonstrated in 
panel (b) of Fig. As expected, in these models the results 
are consistent with observations at an ionization param- 
eter an order of magnitude higher than that required for 
the MF ionizing spectrum. Thus, we can reproduce the 
observed column densities with Sub-Solar metallicity and 
Solar abundance ratios. 
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Fig. 7. Predicted column densities of different species as a function of ionization parameter for log N(H i) (cm -2 ) = 17. 
Results given in panels (a) and (b) are obtained for the MF and BAL spectrum with metallicity 0.1 Zq and 0.5 Zq 
respectively. The results for MF and BAL ionizing spectrum with metallicity 0.1 Zsb of the Starburst enrichment 
(see Table EJl are given in panels (c) and (d) respectively. In each panel the horizontal dotted lines (from bottom to 
top) mark the limits on the column density of Si n , C n , Si iv and C iv obtained from our observations. While the 
former two are lower limits the latter two are upper limits. The verticle dashed lines give the limits on the ionization 
parameters based on the X n/X iv ratios and the observed ratios of N v and Si iv . The verticle shaded region gives 
the ionization parameter range over which the observed column density of N v is reproduced for the chosen nitrogen 
metallicity. 



It is usual in the case of the BLR clouds to explain an 
enhanced nitrogen abundance by rapid star formation nu- 
cleosynthesis. We consider such a chemical enrichment 
in our models using option "starburst abundances" in 
Cloudy (see Table 01 . The details of the models are given 
in Hamann & Ferland (1993). These models have prob- 
lems in reproducing the observations when a MF ionizing 
spectrum is considered (see panel (c) in Fig. 0). However, 



the observed ratios are reproduced if the BAL spectrum is 
used instead (see Tableland panel (d) in Fig.0) for over- 
all metallicity > 0.2 and log N(H i) (cm~ 2 ) in the range 
16 to 17. 

It is interesting to note that for all the models consid- 
ered here the total hydrogen column density is at least 
an order of magnitude less than the amount required (i.e. 
A r (H(total)) > 5 x 10 23 cm" 2 ) to explain the lack of X-rays 
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in this QSO by absorption of a standard QSO spectrum 
(Table|3J). We also run models by fixing the total hydrogen 
column density to be A r (H(total)) > 5 x 10 23 cm -2 . For 
MF ionizing spectrum we do not find any model that can 
reproduce the constraints noted in Tabled Thus, we can 
conclude that either SDSS J160501.21-011220.0 is intrin- 
sically X-ray weak or, if there is any Compton thick screen 
between the quasar and us, it has nothing to do with the 
gas that is associated with the red component. 
The similar characteristics of the blue component together 
with the low 7V(H i) column density of the system strongly 
suggests that the blue component as well is not Compton 
thick. 

4.2. Formation redshift of the blackhole 

We estimate the mass of the black hole from the Eddington 
accretion. For this we estimate the rest frame luminosity 
at 2500 A using the low dispersion spectra obtained by 
the SDSS group (Fan et al. 2000) to be L„=1.37x 10 31 
erg s _1 Hz -1 . Applying bolometric correction as suggested 
by Elvis et al. (1994), we obtain a bolometric luminosity 
of Lfc o ;=9.21xl0 46 erg s _1 . Assuming Eddington accretion 
we derive the black hole mass, Mbh = 7.1 xlO 8 ?? -1 M Q 
(see Haiman & Loeb 2001). Here, rj is the efficiency pa- 
rameter (i.e. accretion luminosity given in the units of 
Eddington luminosity). If e=L{, ;/Ms#c 2 is the radiative 
efficiency for a mass accretion rate Mbh then the natural 
e-folding time scale for the growth of a single seed can be 
written as 

t = M BH /Mbh = 4.0 x lQ 7 e .iT V (3) 

It will therefore take about ln(7.1 x IO^^Mq/IO M ) = 
18.1 e-folding times (or t > OJeo.i?? -1 Gyr) for the black 
hole to grow to the above estimated mass from a stellar- 
mass seed of 10 M©. For the assumed value of e « 0.1 and 
r\ pa 1.0 this corresponds to a formation epoch of the black 
hole close to z ss 11 for the cosmological model considered 
here. 

5. Summary 

We have presented the analysis of broad absorption lines 
(BALs) seen in the spectrum of the z cm ~4.92 QSO SDSS 
J160501. 21-011220.0. Our high spectral resolution UVES 
spectrum shows two well detached absorption line systems 
at z a bs= 4.685 and 4.855. The system at z a bs= 4.855 cov- 
ers the background source completely suggesting that the 
gas is located outside the broad emission line region. On 
the contrary the system at 2 a b s = 4.685, which is redshifted 
on top of the quasar continuum, has a covering factor of 
the order of 0.9. Physical conditions are investigated in the 
BAL system at z a b s = 4.855 using detailed photoionization 
models. The observed H I column density and the limits on 
C II and Si n absorptions suggest that log iV(H i) (cm -2 ) is 
in the range 16-17. The observed column densities of N V , 
Si iv and C iv in the z a bs= 4.855 component require, un- 
like what is derived when analysing broad emission lines, 



that nitrogen is underabundant by more than a factor of 3 
compared to silicon if the gas is photoionized by a typical 
QSO spectrum. Thus, if the gas is ionized by a standard 
MF spectrum, the chemical enrichment of the cloud is dif- 
ferent from that required by emission line clouds. We show 
however that the relative supression in the N v column 
densities can be reproduced for Solar abundance ratios or 
abundance ratios typical of rapid Starburst nucleosynthe- 
sis if we use an ionizing spectrum that is devoid of X-rays. 
Thus if the composition of BAL is like that of the emission 
line regions it is most likely that the cloud sees an ionizing 
spectrum similar to what we observe from this QSO that 
is strongly attenuated in the X-rays. This is consistent 
with the fact that none of our models have high enough 
Compton optical depth to be able to remove X-rays from 
the QSO. Similar arguments lead to the conclusion that 
the system at z a b s = 4.685 as well is not Compton thick. 
Using simple Eddington arguments we show that the mass 
of the central black hole is ~ 8 x 10 s M©. This suggests 
that the accretion onto a seed black hole must have started 
as early as z ~ 11. This gives a typical formation epoch 
for the host galaxy of the QSO. 
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